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Abstract 

Thermal activation is mediated by field configurations that correspond to 
saddle points of the energy functional. The rate of probability flow along the 
unstable functional directions, i.e the activation rate, is usually obtained from 
the imaginary part of a suitable analytic continuation of the equilibrium free 
energy In this note we provide a real-time, non-equilibrium interpretation of 
this imaginary part which is analogous to the real-time interpretation of the 
imaginary part of the one-loop effective potential in theories with symmetry 
breaking. We argue that in situations in which the system is strongly out 
of equilibrium the rate will be time dependent, and illustrate this with an 
example. 
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The analysis of thermal activation has always been an important and interesting topic 
in different areas in physics. Recently, however, this subject has assumed even greater 
importance due to the observation that there are field 

configurations in the standard model (so-called sphalerons[]T], fj) that mediate 
baryon number violating transitions which are unsuppressed at high temperature. Need- 
less to say, this has important implications for the evolution of baryon number in the 
early universeQ. Another reason why there has been a rekindling of interest in the 
evolution of metastable states in the early universe is due to the development of viable 
models of inflation (extended inflation ||) that go back to Guth's0 idea of ending the 
inflationary era via false vacuum decay. In short, there are good reasons for making sure 
that the decay of metastable states at finite temperature is understood from as many 
different perspectives as possible. 

The motivations for this note are twofold. First, we offer a real-time interpretation 
of the standard equilibrium calculation of the activation rate based on the imaginary 
part of an analytically continued free energy. Here we draw an analogy with the case 
of the imaginary part of the effective potential in theories with spontaneous symmetry 
breaking as studied by Weinberg and Wu@. Second, we will make the point that in an 
out of equilibrium situation, different initial states can make for significant changes in the 
activation rate. We show by an explicit example that in an out of equilibrium situation 
the activation rate must be time dependent. 

Two of the seminal works on this topic are those by Langer|| and Affleck [[J. It will be 
instructive to review the concepts and assumptions involved in these calculations, since 
our results will provide a different interpretation of these standard results, as well as an 
extension of them. 

Langer|8], [10] develops a Fokker-Planck type equation for the probability of finding 



the system in a given configuration at time t. There are two terms in this Fokker-Planck 



equation. The first term is deterministic in nature, corresponding to the Liouville time 
evolution in phase space. The other is a stochastic term and represents the diffusion or 
dissipation arising from the coupling to a heat bath. The probability density in phase 
space obeys a continuity equation and the associated current gives the flow of probability. 
In a metastable situation, there is a particular configuration that corresponds to a saddle 
point of the energy functional, and probability flows along the unstable direction from 
the metastable phase towards the stable phase. The activation (decay) rate is obtained 
as the total current passing over the saddle point along the unstable direction divided by 



the initial population of the metastable state [|1CJ . 

There are two very important assumptions in Langer's original work. The first is 
that a steady state solution is reached in which a source in the metastable state contin- 
ually feeds probability (or particles) and a sink in the stable phase removes probability 
at a steady rate. Secondly, the metastable state is supposed to be in local thermody- 
namic equilibrium. This assumption implies that there is a wide separation between the 
equilibration (relaxation) time scale and the time it takes for the state to lose a sig- 
nificant amount of probability. This in turn translates into the statement that there is 
fairly strong dissipation JlO|. Under these assumptions, there is no memory of the initial 
preparation of the metastable state. 

Similar assumptions are implicit in Affleck's calculation ||. The rate is defined as the 
Boltzmann average of the decay rates of the quantum states that, for energies less than 
the barrier height, are standing waves in the metastable well. For energies higher than 
the barrier height, these are waves incident from the left, reflected and transmitted at 
the barrier. The rate is then calculated as: 



POO 

r = Zq 1 / dEp(E)T(E) exp(-/3 E), 
Jo 
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where p(E) is the density of states at energy E, T(E) the decay rate for states with 
energy E, and (3 = 1/T. The rate is normalized using the partition function Zq of a 
harmonic oscillator centered at the metastable state. As in Langer's calculation, this 
corresponds to a steady flow of probability across the barrier. 

Both of these calculations (and those that have built on them) assume that the 
state under consideration is one in which probability is being fed continously into the 
metastable well in order to replenish the probability that is flowing out and over the 
activation barrier, thus ensuring a steady state. Alternatively, the assumption is that 
the system is in quasi-equilibrium and "slowly" leaking probability on times scales much 
longer than typical scales of the system. However, this implies that the system is looked 
at during times for which all transient effects arising from the initial preparation have died 
out and that the probability inside the metastable well has not been depleted appreciably. 



In particular the underlying assumption! 1C] is that the equilibration time inside the 



metastable well is much shorter than the inverse rate of depletion of the population in 
the well. 

These authors showed that under these (implicit) assumptions, the rate obtained is 
the same as that obtained from the imaginary part of the analytically continued free 
energy: 

r = ImF (2) 

where Q is the unstable frequency at the top of the barrier. 

This description offers a fairly reliable method of calculation in a wide variety of 



experimental situations in macroscopic systems^^, in which there is a strong coupling 
to an environment (heat bath) and only a few degrees of freedom are relevant. 

However, the setting we have in mind is that of a weakly coupled quantum field theory, 
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in which the separation between system and bath is somewhat blurred. In particular, 
at the onset of a first order phase transition, several time scales are relevant such as the 
equilibration rate, expansion rate (in a cosmological situation) and the activation rate. 

Given that this calculation of the rate involves assumptions which in any typical 
quantum field theory will be extremely hard to justify at the microscopic level without 
invoking some phenomenological approximation, we want to investigate an alternative 
real-time description of the process of thermal activation. We also envisage situations in 
which the system may be strongly out of equilibrium, and in which some of the above 
assumptions are not likely to apply. 

To illustrate our discussion with a definite example, consider the case of a scalar 
field theory described by a potential with a local metastable minimum and a global, 
stable minimum. In this case, there is a static spherically symmetric solution to the 
classical field equations that resembles a droplet of the stable phase inmersed in a sea 
of the metastable phase[|S], [T7]]. This solution is characterized by the following collective 
coordinates: translational ("zero modes") corresponding to the position of the droplet 
and another that represents the radius of this droplet (although technically speaking this 
latter does not correspond to a "zero mode"). The static solution to the field equation 
corresponds to a particular value of the radius of the droplet, the critical radius, at which 
the surface energy is balanced by the volume energy. 

Fluctuations around this particular solution are characterized by zero frequency modes 
(corresponding to translational invariance), a negative eigenvalue corresponding to small 
fluctuations of the radius of the droplet, and a spectrum of positive frequencies. Thus 
this configuration corresponds to a saddle point in functional space. The unstable mode 
describes the instability towards collapse or growth of the droplet for radii smaller or 
larger than the critical radius, respectively, when either the surface energy or the volume 
energy start to dominate. 
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Quantization around the classical solution corresponds to treating the position of the 
droplet as a collective coordinate and expanding |p"T|, |T2| , |T5| |H| , p75f| 



$(f, t) = $ d (f - r (t), = R*) + J] Qi(t)Mf-f (t), R = R*), (3) 

where f D is the position of the droplet, R is the radius and R* is the critical radius 
for which the energy functional has an extremum. The mode functions fi(f) are the 
eigenf unctions of the fluctuation operator around the static droplet configuration and are 
chosen to be orthogonal to the "zero modes" [HJ. There is a coordinate Q u associated 
with an unstable mode with mode function f u (r) ~ d^ d /dR\ R=R * and negative frequency 
— Q 2 . To lowest order in the semiclassical expansion, the Hamiltonian operator becomes 



fi 2M 2 2 V 



p2 1 
V I ,2^2 



(4) 



2 2 

with being the energy of the critical droplet, M is given by the normalization of 
the zero modes and ujv are the stable frequencies. To this order in the semiclassical 
approximation, P is the total momentum of the system. When the radius of the droplet 
is much larger than the correlation length (thin wall approximation) 

Att 

E(R) ss AitR 2 a R 3 AS (5) 

3 

with a being the surface tension given by the gradient terms and AS the energy (or free 
energy) difference between the globally stable and the metastable lower energy (or free 
energy) states for the homogeneous configuration. In this case the radius of the critical 
droplet is given by 

V « % (6) 

In a saddle point approximation of the equilibrium partition function, only two ex- 
tremum configurations are usually considered. The first corresponds to the metastable 



minimimun, while the second corresponds to the droplet configuration (saddle point). 
The partition function is thus written as the sum of these two independent terms. 

Formally, the partition function evaluated at the saddle point corresponding to the 
droplet does not exist because of the presence of the unstable (negative frequency squared) 
oscillator. However, it may be defined via an analytic continuation, in which case it at- 
tains an imaginary part which in three spatial dimensions is given by: 



±iV 



MT~\ 3 / 2 exp[-M 



2tt 



2 sin(ft/2T) A J# A 2 sinh(^/2T) ' 



(7) 



where the ± arises from the direction of the analytic continuation. The volume factor V 
and the T 3 / 2 factor arise from the "free particle" zero mode, the 1/ sin(fi/2T) from the 
unstable mode and the product over /' is only for the stable modes. Writing the partition 
function Z as: 



1 + 



Z G exp 



with Z a the partition function of independent harmonic oscillators at the metastable 
minimum of the potential, choosing the negative sign for the analytic continuation (so 
as to obtain decaying exponentials), and using the identification (Q) finally leads to the 
usual formula for the rate per unit volume: 



r 

V 



MT 



exp[ 



n,2sinh(^/2T) 



(9) 



2vrsin(fi/2T) Ui> 2 sinh(^//2T) ' 
Here lo% are the frequencies at the metastable minimum and iov are the stable frequencies 
in the perpendicular directions at the saddle point. 

In functional space, the direction of steepest descent at the saddle point corresponds 
(locally) to the unstable coordinate Q u . For the purposes of analyzing the decay of 
a metastable state and activation over the barrier (near the top of the barrier), only 



7 



this quantum mechanical coordinate is important. Thus along this unstable direction, 
the problem becomes effectively that of a quantum mechanical inverted oscillator near 
the top of the barrier. This approximation implies that the temperature is such that 
ti)t<T< E(R*), where u k are the typical frequencies at the metastable well. 

Another situation in which the free energy (or ground state energy) acquires an imag- 
inary part is in the calculation of the one-loop effective potential for homogeneous config- 
urations either at zero or finite temperature in the case when the tree level potential V(<&) 
has a double well structure 0. In this case the fluctuations of wave vectors k 2 < V"(<&) 
are unstable since they see an inverted harmonic oscillator, just like the unstable coor- 
dinate above. These modes are treated via an analytic continuation similar to the one 
performed above, giving rise to an imaginary part for the one loop effective potential (or 



free energy) for homogeneous configurations [1C 



Weinberg and Wu0 argued that this imaginary part of the effective potential can 
be interpreted as the decay rate of a gaussian state peaked at the "top of the inverted 
oscillator" for these unstable modes. The rate of spread of this gaussian wave packet was 
identified with the imaginary part of the one-loop effective potential. 

Further analysis JL6] shows that this rate also determines the rate of growth of fluc- 



tuations which is determined by the width of the gaussian state initially prepared to be 
centered at the top of the inverted parabola at the maximum of the tree level potential. 
This analysis was also carried out at finite temperature, in which case the analytically 
continued free energy has an imaginary part[l(J. These unstable modes are responsible for 



the early stage dynamics of a second order phase transition and for the non-equilibrium 
evolution during the phase transition. The regions in which the field becomes correlated 
grow in time as the system becomes more and more correlated. Initially the unstable 
modes grow exponentially but eventually the non-linearities set in, slowing down the rate 
of growth. Finally, at a later time the growth stops. The system is now described as 
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having large correlated regions in which the field inside is close to one of the minima of 
the (effective) potential. This is the mechanism of "spinodal decomposition" |TjJ. Thus 
the imaginary part of the one-loop effective potential (or alternatively at finite temper- 
ature the free energy) disguises a time dependent situation. Trying to study this time 
dependent situation via an equilibrium description leads to treating the unstable modes 
via an analytic continuation and an imaginary part. 

The physical difference between the mechanism of spinodal decomposition and ther- 
mal activation is that in the former there are no free energy barriers, and small amplitude 
long wavelength fluctuations become unstable and grow. By contrast, in thermal activa- 
tion, there is a free energy barrier to be overcome. Thus, a large amplitude configuration 
must be present in the bath with a radius that is bigger than the critical radius. Once 
this configuration is created, however, it is the instability towards growth of this bubble 
that drives the phase transition. The probability of finding these configurations in the 
bath is determined by the Arrhenius- Kramer exponential factor explicit in Z an above. 
When a critical bubble is formed, the system is found at the top of the barrier, that 
is at the saddle point. In the gaussian approximation, the situation corresponds to a 
gaussian density matrix (the thermal equivalent of the gaussian wave packet studied by 
Weinberg and Wu) centered at the top of the barrier, and for the unstable coordinate 
corresponds to an inverted harmonic oscillator. From this point onwards, the analysis 
of the evolution, is then similar to that of Weinberg and WuJJJ and Boyanovsky and de 
Vega||16||, with the difference that the initial state is a thermal density matrix rather than 



a wavepacket. 

We see then that the imaginary parts of the free energy and/or the effective potential 
usually found in the calculation of the decay rates are a signal of non-equilibrium, real- 
time evolution. 

Thus motivated by the situation with the one-loop effective potential we now present 



a real-time evolution of an initial state that corresponds to a one- dimensional quan- 
tum mechanical system describing the unstable coordinate Q u . The reduction of the 
multidimensional problem to one collective coordinate is justified because this unstable 
coordinate determines (in a neighborhood of the saddle point) the direction of steepest 
descent in functional space. This is the direction along which probability will flow. 

The procedure is straightforward. Start with an initial density matrix p(t = 0) and 
then evolve it in time via either the Liouville equation: 



ih^- = [H,p(t)\, (10) 



or via the solution to this equation: 



p(t) = exp(-~Ht)p{0) exp(ifft). (11) 
Here H is the Hamiltonian for the unstable coordinate 

h u = y~ \ n2 Ql ( 12 ) 

Given the density matrix as a function of time, we can look at its position space 
representation p(Q u , Q' u ',t) = (Q u \p(t)\Q' u ) . The current along the unstable direction is 
then found via: 

J(Qu,t) = |(^- - ^r)p(Qu,Q' u ,t)\ Qu=Q ,. (13) 

Evaluating this current at the saddle point will then give us the transition rate (or 
activation rate) over the barrier. 

Generalizing the analysis of Weinberg and Wu|7] and Boyanovsky et. al. ||16|| , we will 



assume that the initial density matrix p(t = 0) (only for the unstable coordinate) is that 
of an upright harmonic oscillator of frequency u) in thermal equilibrium at temperature 
T = 1/(3 given by|TJ[: 
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p(Qu,Q' u ;t = 0) =^(0)exp 



2hsmh((3?uv) 

W) - 



(Ql + Q' 2 ) cosh(fihu) - 2Q U Q' U ] (14) 



(15) 



/ 27rh sinh((3huj) 

The coordinate space expression for the time evolved density matrix is, for t > 0: 

/% % 
dy dy'(Q u \exp(-~H u t)\y)p(y,y';t= 0)(y'\exp(-H u t)\Q' u ) (16) 

The propagators (x\ exp(±~H u t)\y) are easy to evaluate by analytically continuing 
the propagator for a standard harmonic oscillator with real frequency [[HJ: 



(x\exp(±-H q t)\y) 



N(t) 



N(t)exp(± 



n 



2frsinh(f2t) 

U/2 



[x 2 + y 2 ) cosh(fit) - 2xy ) (17) 



2rrih sinh(fit) 

It can be easily checked that these propagators are solutions to the evolution equation 
with the proper boundary conditions. We can now compute the density matrix as a 
function of time, as well as the current. We have verified that the resulting density 
matrix is a solution of the Liouville equation ([LT]) with the initial boundary condition 
given by (|T4"D, thus confirming that the analytically continued propagators give the correct 
answer. 

Rather than write down the density matrix, we consider the probability density 

p{Qu,t) = p(Q u ,Q u ;t): 



p{Qu,t) 

a(t) 



Af(t) exp(-Q 2 u /2a(tf 



h 



2iutanh((3huj/2) 
a(0)- 



cosh 2 (Ot) + ^-sinh 2 (^t) 



Af(0) 



a{t) 
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Thus, at long times, the probability to find the density matrix of the unstable coor- 
dinate centered at the saddle point is 



p(0,t) 



— exp[— Qt] 



(19) 



[1 +LU 2 /tt 2 }- 

At zero temperature this is similar to the result obtained by Weinberg and Wu[[/J for a 
particular unstable mode by identifying Q with the unstable frequency for that particular 
mode. 

The probability current (|T3| ) evaluated at the saddle is zero because the density matrix 
spreads symmetrically around the saddle point along the unstable direction. However, 
we may still define a rate as 



dln\p(0,t)] 
dl 



(20) 



This expression will be valid for t ^> Q^ 1 . Integrating this rate in all functional directions 
perpendicular to Q u to obtain the total probability at the saddle point and dividing by the 
partition function for the harmonic oscillators at the metastable well, we finally obtain 
an alternative definition for the rate per unit volume: 



r 

-± = n 
v 



MT 
~2^~ 



3/2 r E(R* 
exp[ — 



n fc 2sinhK/2T) 
III' 2 sinh(c^/2T) 



(21) 



which differs from ([]) by the prefactor. 

Each contribution to the rate above has a very simple and clear real-time interpreta- 
tion. The situation considered corresponds to the system being described as a gaussian 
density matrix centered at the saddle point. Along the unstable direction the density ma- 
trix (similarly to a wave packet) spreads at a rate determined by the unstable frequency 
Q. Thus the rate (|2lf) corresponds to the total rate of spread of the initial density matrix 
(integrated over all the perpendicular directions) divided by the total probability in the 
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metastable well. The underlying assumption here is that the system has reached the 
saddle point and remains there in quasi-equilibrium and that the depletion of probabil- 
ity at the saddle is a consequence of the spreading of the density matrix, although the 
statistical average of the collective coordinate (radius of the droplet) remains at the top 
of the barrier. 

We may draw the following conclusion from this analysis. The calculation of the rate 
via the imaginary part of the free energy, which in Langer's work is derived from a steady 
state assumption, corresponds to the assumption that the statistical density matrix is 
centered at the saddle point, and that the expectation value of the radius of droplet is 
the critical radius. Hence, just as in the case of the one-loop effective potential, the 
imaginary part of the free energy is disguising an intrinsically time dependent situation 
and it will determine the rate of growth of fluctuations (width of the gaussian density 
matrix) 0, [16|]. This can be easily understood from the fact that the fluctuation in the 
unstable coordinate is given by 



(Ql(t)) oc a\t) (22) 

We can apply the same analysis for a very different situation, that of the system initially 
away from the saddle and for which a steady state assumption is not applicable. Thus 
we now consider activation as an initial condition problem and look at the situation in 
which the average radius of droplets is smaller than the critical. This corresponds to 
an initial situation far from equilibrium and not in the steady state. Such a condition 
corresponds to the situation studied in a simplified 1 + 1 dimensional field theory model 



by Boyanovsky and de Carvalho |2C 



Again we will only concentrate on the unstable coordinate Q u , and consider the case 
in which the initial density matrix for this collective coordinate is that of an upright 
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harmonic oscillator of frequency u that is displaced a distance Q uo < from the origin 
at the initial time t — and with zero average of the momentum conjugate to this 
coordinate: 



p(Qu,Q' u ,t = 0) = jV(0)exp 
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2h smh(phu 

2(Q u -Quo)(Qu-Quo)}} 



((Ql - Qlo) + (Qu - Qlo)) cosh(phu) 

(23) 



with the same normalization factor as in (jl4j) . 

We now solve the Liouville equation along the same steps used in the previous 
case using the real-time propagators. We obtain the probability as a function of time 

p(Qu,t) = p(Q u ,Q u ;t) as: 



p(Q u , t) = Af(t) exp -{Q u - Q uo cosh(^)) 2 /2a(t) ; 



(24) 



with a{t) and N{t) given by (|i~8p. The current (|T3| ) along the unstable direction evaluated 
at the saddle point {Q u = 0) is found to be 



■UQu = 0;t) = — Af(0)\Q uo \A(t)exp [-taah(Phu/2)B{t)] 



(25) 



with 



A(t) 
B(t) 



sinh(fit) 



cosh 2 (fit) + to 2 /n 2 sinh 2 (fit) 



3/2 



luQ 



2 

uo 



to 



1 + '—tanh 2 (nt) 



(26) 
(27) 



The current is shown in figure 1. The probability evaluated at the saddle point is 



p(Q u = 0,t) =Af(t)exp [-(Q uo cosh(nt)) 2 /2a(t) 2 
with a(t) and Af(t) given by (|TJ|). At times t > Q' 1 , the rate 



(28) 
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coincides with (|20| ) and only reflects the spread of the initial density matrix. One can use 
this rate as an alternative definition in real time, leading again to the result ([21]) after 
integrating along all the perpendicular directions in functional space. 

The depletion of probability at the saddle clearly reflects two physical effects. The first 
is the spread of the density matrix determined by cr(£), while the second is the "rolling" 
of the unstable coordinate down the inverted harmonic oscillator potential. The same 
features give the total current at the saddle. However, whereas the coordinate rolls down 
with a time dependence Q uo cosh(f2t), the spread of the density matrix increases much 
faster and is given by o~(t). The two effects thus combine to give a positive current along 
the saddle point. We see then that despite the fact that bubbles with radii smaller than 
the critical radius collapse, they still contribute a positive current over the barrier. The 
total current along the unstable direction, evaluated at the saddle point is obtained by 
integrating the current ( p5|) along all the perpendicular directions ||, and restoring the 
exponential Arrhenius factor. It is normalized by dividing by the partition function of 
the metastable well. This procedure then yields: 



<Jt,u — Ju{Qu — 0;£) 



MT 



3/2 



E(R*) n fc 2sinh(^/2r) 
P[ T J riz'2sinh(^/2T) [ ^ } 



2tt 

where the factors were obtained as in ( ^IP and J U (Q U — 0; t) is given by (|25j) . 

There are two questions that require an answer at this point: i) which one is a suitable 
definition of the rate, ii) what is the range of validity of the approximations involved in 
obtaining the results presented above. We answer these in turn: 



i) We have seen that the definition of the rate given by (|2T| ) coincides (up to the 
prefactor) with the oft quoted result in the literature given by equation (Q). This defi- 

15 



nition corresponds to describing the process of activation by assuming that the system 
is described by a gaussian density matrix in terms of the small fluctuations around the 
critical droplet, and that the activation rate is given by the rate of spread of the probabil- 
ity for a critical droplet integrated along all perpendicular directions in functional space. 
Although this seems a sensible definition of the rate, it is not the definition of the rate of 
change of probability for finding the system localized in the metastable state. This latter 
definition corresponds to calculating the total probability current at the saddle point as is 
obtained either from quantum mechanics or via the Fokker-Planck description P, |l(J and 
corresponds to the current ( p9|) calculated above. The steady state assumption implicit 
in Langer's treatment thus leads to the identification of the two. However, we expect 
that in a strongly out of equilibrium situation, for which a steady state assumption is 
not valid, or during times for which the details of the initial state are important, the 
correct definition of the rate, that is, the probability current at the saddle, will yield a 
time dependent rate which will be sensitive to the initial conditions. In the case we have 
described above, this corresponds to sensitivity to the initial value of the coordinate Q u . 
ii) There are several approximations invoked in this (and most other) analysis. First, the 
quadratic approximation around the critical droplet configuration. This approximation 
will be justified insofar as the higher order terms in the expansion (cubic and quartic 
in terms of the coordinates Qi) can be neglected. Given that the imaginary part of the 
free energy describes a time dependent process in which the width of the probability 
distribution grows and therefore the fluctuations are growing, this approximation will 
break down at some time that will be given by the details of the particular problem, 
such as couplings and structure factors of the stable mode functions fi{r). In order to 
identify the rate fl2T|) we had to assume that t ^> fl -1 , for small Q (as is expected in the 
thin- wall approximation). This approximation may imply times longer than the regime 
of validity of the gaussian approximation, in which case the rate may depend on time. 
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A quantitative analysis of the validity of the approximations will require a more detailed 
knowledge of the potentials, mode functions couplings etc. and will have to be done for 
each particular problem individually. 

We started this work by using examples such as first order inflationary models and 
baryon number violation via sphaleron mediated decays to motivate the discussion on 
thermal activation. It is our impression that there is certain amount of discomfort and 
uneasy trust on the usual calculation of the activation rate via the imaginary part of an 
analytically continued equilibrium free energy. We have drawn the analogy between the 
case of thermal activation and the case of the one-loop effective potential for theories 
with spontaneous symmetry breaking at tree level. Both calculations involve gaussian 
approximations around unstable configurations (inverted harmonic oscillators) and both 
calculations obtain an imaginary part by analytically continuing the result for upright 
harmonic oscillators. 

In both cases, the imaginary part of the free energy is associated with the growth of 
the width of an initially prepared gaussian state (density matrix) and signals the growth 
of fluctuations as a consequence of the instabilities. 

We have shown that a real-time description of this process leads to a definition of the 
rate that, up to prefactors, is similar to that obtained via the recipe using the imaginary 
part of the free energy. 

We also argued that in general, non-steady state or strongly out of equilibrium situa- 
tions the correct definition of the rate, as the total probability current flowing across the 
saddle point along the unstable direction will lead to a rather different result which may 
in general be time dependent. We illustrated this situation with a non-equilibrium case 
with an initial situation in which the mean value of the collective coordinate represent- 
ing the radius of the droplet has an expectation value which is smaller than the critical 
radius. Time dependent rates are common in activated events in condensed phases pTf 
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where transient behavior after an initial preparation is observed (see also [pP| ) . 

Although we cannot make general claims as to what is the correct rate in a given 
theory, our goal here was somewhat more restricted: to provide a dynamical, real-time 
interpretation of the imaginary part of the free energy and at the same time investigate 
scenarios for which a calculation of the rate based on this analytic continuation may not 
yield reliable results. 

We are currently involved in calculations of decay rates in theories involving sphalerons 



using the real-time formalism developed here|22[ as an initial condition problem. Specif- 
ically, if the phase transition occurs via a quench in such a way that the system does not 
have time to equilibrate around the minima, thermal activation will occur strongly out of 
equilibrium and dependence on the initial state and time dependent rates are expected. 

The issue of thermalization in the wells requires a deeper understanding of the time 
scales involved, in particular in inflationary scenarios. 

While we do not yet have all the answers we need to fully understand what changes 
the time dependence of the rate will bring, we may speculate. We have considered 
thermal activation here rather than under the barrier tunnelling. This makes the range 
of applicability of our calculation to inflationary models somewhat suspect. The reason 
for this is that once inflation sets in, the temperature of the heat bath will decrease 
rapidly, turning the problem into a zero temperature one. Even in this case, however, we 
should expect the rate to be time dependent. The basic change in our calculation is that 
the paths used to compute the propagators required to evolve the density matrix in time 



will be different. Essentially, one must do a WKB approximation of the propagators p3 

The situation in extended inflationary models is somewhat trickier to assess, since 
the nucleation rate is already time dependent in most of these models due to the time 
evolution of the Jordan-Brans-Dicke field in them 123 



In the case of the sphaleron, the question of whether the sphaleron interactions are 
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in thermal equilibrium (which is crucial in terms of determining whether a baryon asym- 
metry can be generated by these interactions), becomes more difficult to assess due to 
the time dependence of the rate. One could imagine that the rate of these interactions 
decreased sufficiently quickly so as to allow them to drop out of local thermal equilibrium 
thus allowing a net B asymmetry to be generated. However, the answer to these and 
other questions will only be found once a deeper understanding of the real-time dynamics 
of the process of thermal activation is understood in these systems. 
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Figure Captions 

Figure 1: The thermal activation rate T(t) as a function of t at fixed temperature. 
Time is measured in units of Q^ 1 . 
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